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The Sample 

The Hg 1 _,Cd Te alloy with x = .1 was prepared by reacting 
the constituent elements in a sealed, f used-si 1 ica, 5-mm i.d. and 
10-mm o.d. ampoule. 

The ampoule was cleaned with aqua regia, distilled water, 
acetone and then methanol , and was baked at 1000 C in a vacuum of 
500 torr (the best available, better than no vacuum). 

The ampoule was then loaded with 1.7835 g Hg first, then 
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1.234 g Te, and finally .108 g Cd , then evacuated to about 10 “ 
torr and sealed-off. 

The sealed ampoule was put in the rocki ng-f urnace. The 
furnace was heated up from room temperature to 760 C (the actual 
temperature is about 90 C higher than this reading from the 
programmable controller) in 10 hours. the furnace was rocking 
while the temperature was held at 760 C for 17 hours, then the 
power and rocking motor were turned off. The ampoule was pulled 
out of the furnace at 520 C. 


The Antimony Reference 

The ampoule used for the reference was identical to that 

used for the sample. The ampoule was cleaned with aqua regia, 

distilled water, acetone, and then methanol, and then baked from 

room temperature to 1050 C at 10 7 torr. 

The ampoule was then loaded with 2.6552 g Sb (same volume as 
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the HgCdTe sample) , then evacuated to about 10 torr. The 
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antimony was melted with a torch, then the ampoule was sealed- 


off. Both the sample and reference ampoules were shortened a 

bit by Robbie to reduce the free volume. 


Temperature Profile 

The Marshal 1 platinum-wound tubular 
vertically. A quartz tube of 30-mm o.d. and 

with fiberfrax and is placed inside a metal 

Controller . 

figure. . 
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furnace is mounted 
27-mm i.d. is filled 
tube as shown in the 


First the temperature profiles at setpoints 650, 900 and 
1100 C were measured with a type-K thermocouple. As shown in 
the plots. x is measured from the bottom of the furnace, and the 
furnace is 24 inches (61 cm) long. Since the peaks of the 
profiles were not consistent, the thermocouple was replaced by a 
type-R thermocouple. Although the profiles were consistent, the 
peak was off center. It was displaced toward the top of the 
furnace by about 3 cm. The flat zone is about 3 cm at 1100 C. 


At lower set points the flat zone broadens slightly, being 4 to 5 
cm at 900 C and about 7 cm at 650 C. 

While the sample and reference were being prepared, the 
temperature profiles described below were measured at a setpoint 
of 900 C. The intention was to find a way to broaden the flat 
zone. A schematic diagram of the furnace and its taps is shown below. 
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a) Shunt between external binding posts 2 and 4 with a resistance 
11, 7.4 and 3.5 ohm respectively. As shown in the plots, putting 
the shunt between taps 2 and 4, makes the profile more peaked, 
and reduces the breadth of the flat zone. 

b> Shunt between taps 3 and 4 with a resistance 3.5, 0.9, 0.4 and 
0 (shorted) ohm respecti vely. From the plots, it shows that 
shunting the center zone (taps 3 and 4) had no effect upon the 
breath of the flat zone except for shifting it to the geometrical 
center of the furnace. The flat zone is essentially the same 
with and without a shunt. 
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c) Several attempts were made to broaden the flat zone by 
shunting taps 4 & 5 and 2 S< 3, in addition to (b) , but these 
efforts resulted in asymmetric profiles which were quite 
unexpected. It is possible that these complications result from 
convection in the vertical orientation or from some peculiarity 
in the windings. 
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CONVECTION IN A VERTICAL AMPOULE DUE TO 
NON-AXIS YMMETRIC TEMPERATURE DISTRIBUTION 

J. Iwan D. Alexander, Jalil Ouazzani and Franz Rosenberger 
Center for Microgravity and Materials Research 
University of Alabama in Huntsville 


1. Introduction 

Convection and segregation in directional solidification and crystal growth by the 
Bridgman-Stockbarger technique have been traditionally treated by assuming axisymmetric 
thermal conditions on the ampoule wall. For most recent work and earlier references see [1], 
It is, however, difficult to achieve such a condition in an experimental setup. Conventional 
heating devices (furnaces) typically have several degrees asymmetry in the temperature 
distribution of a plane normal to the geometrical symmetry axis. This can be significantly 
improved upon using heat pipes. But the Bridgman technique requires axial temperature 
gradients, i.e. transition zones between two quasi-isothermal zones (heat pipes) at different 
temperatures. Furthermore the temperatures required for the melt growth of most 
technologically important materials are high enough that radiative heat transfer to the ampoule is 
significant, if not dominant. Then, any axial misalignment of the ampoule, even in a perfectly 
axisymmetric gradient, will lead to a non- axisymmetric temperature distribution on the ampoule 
wall. This is due to the different view factors of the ampoule surface elements in a plane 
normal to the ampoule's axis. 

Any deviation from an axisymmetric temperature field on the wall of a vertical ampoule 
represents a horizontal temperature gradient. The horizontal density gradient that results 
from thermal expansion in the melt under this condition must lead on Earth to some 
bouyancy-driven convection, no matter what the axial (vertical) temperature distribution that 
is imposed on the melt. It is the purpose of this study to estimate the magnitude of such 
convective flows for conditions representative of the MSFC mercury-cadmium-telluride (MCT) 
Bridgman setup. 
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2. Model and Method of Solution 

The model consists of a closed vertical cylinder of radius R=0.25cm and height 
H=0.5cm. A vertical, convectively stabilizing temperature difference of 4.4 *C is imposed on 

the cylinder with the temperatures of the lower and upper plates, respectively, being kept 
constant at ^=958.1 and T 2 =962.5 *C. 



TABLE 1 


Location Temperature [°C] 


A 

958.115 

B 

958.346 

D 

958.807 

E 

959.038 

F 

959.268 

G 

959.499 

H 

959.729 

I 

959.959 

J 

960.190 

K 

960.420 

L 

960.651 

M 

960.882 

N 

961.112 

0 

961.343 

P 

961.573 

Q 

961.804 

R 

962.034 

S 

962.265 

T 

962.497 


y 

(b) 


Fig. 1. Defining sketch, (a) Coordinate system and wall temperature (isotherms) distribution, 
(b) Table of temperature values for isotherms in (a). 


On the side wall, as depicted in Fig. 1, a temperature distribution of the form 


zR 


n z 


T wal i= T x + (T 2 - T x ) 1-cos cos 9 


2H 


is imposed. This amounts to a maximum horizontal temperature difference of 2.2 °C. The 
cylinder is completely filled with a fluid of Prandtl number Pr=10‘ 2 (or thermal diffusivity k = 
OxlO' 1 cm 2 /sec and kinematic viscosity v=1.3xl0' 3 cm 2 /sec) and thermal expansion 
coefficient (3 = 2.45x10^ K' 1 . These parameter values are representative of molten germanium 
[1]; corresponding values for MCT alloys are, unfortunately, not known with sufficient 
accuracy at this point The acceleration of gravity is assumed to act downward, parrallel to the 
cylinder axis. For simplicity, its value is assumed to be g=10 3 cm/sec 2 . 


It is assumed that the melt behaves as a Boussineq fluid, 
governing the transport of momentum and heat are then 


The basic equations 


— + (grad u)u = - — gradp + vAu - gk, 

Po Po 


div u = 0, 


( 1 ) 

( 2 ) 


— +u»gradT =kAT, (3) 

3t 

where u, p, and T are the fluid velocity, reduced pressure and temperature, v is the kinematic 
viscosity, p G is the (constant) reference density, k is the unit vector in z-direction and Ap/p 0 is 

the change in density due to temperature. On the boundaries we impose, in addition to the 
above thermal conditions, the no-slip condition u = 0. 

The system of equations (1-3) was solved for these boundary conditions employing the 
finite element code FIDAP [2,3]. The computational domain was discretized with the mesh of 
8 node bricks depicted in Fig. 2. The non-isometric grid in physical space was transformed to 
a grid of isoparametric bricks in logic space. For each of these elements the velocity 
components and the temperature were approximated using trilinear interpolation functions. The 
pressure was approximated using a piecewise constant discontinuous pressure approximation 
(penalty method)[4]. For the steady-state algorithm we used a Newton-Raphson method 
combined with Quasi-Newton updating [5,6]. 
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Fig. 2. Mesh system used for finite element 
solution. 


3. Results 

The major result of our modelling is that the imposed temperature distribution leads to 
convective flow of significant strength. As illustrated by the velocity field in the symmetry 
plane (y=0) presented in Fig. 3a, one convective roll is established. Note that the maximum 
velocity magnitude is in excess of 0.25 cm/sec. Figure 3b shows the resulting temperature 
distribution in the symmetry plane. The flow has strong three-dimensional features, as is 
illustrated by the velocity plots for various cross-sectional planes that are presented in Figs. 4-6 
together with the corresponding temperature distributions. Note that the velocities in all these 
cross-sections are of comparable magnitude. 
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Fig.5. Plane x=0. Velocities at points in 
plane, (a) Vector representation, (b) 
velocity surface. Maximum vector 
plotted corresponds to 0.23 cm/sec. 
(c) Isotherms in plane; for values see 
Table 1. 
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Fig.6. Midplane z = H/2. Velocities at points in 
plane; (a) vector representation, (b) 

velocity surface [rotated -45° with 
respect to (a)]. Maximum vector plotted 
corresponds to 0.27 cm/sec. 

(c) Isotherms; for values see Table 1. 


4. Discussion and Conclusions 

In view of the astounding magnitude of the convective velocities obtained, we have 
used the basic solution of Hart's sideways diffusion model (setting the solutal gradient to zero) 
to calculate the maximum velocity associated with a horizontal temperature difference 
maintained between two vertical planes [7]. One would expect that Hart's model, due to less 
flow impedance by no-slip surfaces, results in a somewhat larger convection velocity than our 
totally confined geometry. This expectation was, indeed, confirmed by the value of 0.83 
cm/sec that we obtained for the maximum velocity, as compared to 0.27 cm/sec from our 
cylindrical model for otherwise similar conditions. 

It should be emphasized that such high convection velocities that result from "small" 
deviations from axisymmetric thermal boundary conditions, exceed diffusion velocities 
characteristic of semiconductor melts by some three orders of magnitude. These convection 
velocities also exceed those obtained in earlier Bridgman modelling work (Grenoble furnace) 
[1], assuming axisymmetry, by one to two orders of magnitude. In addition, the convective 
flow resulting from the asymmetry condition extends throughout the melt. This is in contrast 
to the highly localized, interfacial flows obtained in the earlier modelling. These localized 
flows arise from the mismatch of the thermal conductivities of melt, crystal and container. 

Our preliminary study needs to be expanded considerably. For instance, the stabilizing 
effect of larger vertical temperature gradients needs to be explored. Furthermore, the effect of 
reducing the severity of the asymmetric thermal boundary condition needs to be investigated by 
imposing flux conditions rather than a fixed temperature profile. But is it probably safe to say, 
even at this early point, that our model shows that symmetric transport conditions in the 
Bridgman technique will be very difficult to realize in the laboratory practice. Future modelling 
must recognize this. 
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